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METHOD FOR FINDING QUOTIENT IN A DIGITAL SYSTEM 


BACKGROUND OF THE INVENTION 


The present invention relates to a method for finding a 
quotient, especially to a method for finding a quotient in a 
10 digital system by signed-digit operation. 

Inherently, division operation is a sequential 
operation* The quotient digits are produced only after the 
sign of the remainder have been detected. As a result, 

15 division operation is much slower than multiplication 
operation. Efforts have been put in speeding up the division 
operation. It is noted that the SRT algorithm (C. V. 
Freiman, "Statistical Analysis of Certain Binary division 
algorithms," Proc. IRE, Vol. 49, Jan. 1961, pp. 91-103; K. 

20 Hwang, Computer Arithmetic: Principles, Architectures, and 
Design, 1979, pp. 222-223) eliminates the restoring 
operations of the partial remainders. Another algorithm 
disclosed by K. Hwang in the aforementioned article confines 
the quotient digits either to be 1 or -1, depending on the 

25 signs of remainders. However, the bottleneck of those 
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algorithms lies in sign detection of the remainder. Fast 
addition algorithms such as CLA (carry-lookahead addition) 
shorten the operation time , but results in complex hardware 
structures. The aforementioned articles of C. V. Freiman and 
5 K. Hwang are hereby incorporated herein by reference. 

Recently, division algorithms based on SD (signed- 
digit) number representation was proposed which is much 
faster than the previous algorithm (S. Kuninobu et al., 

10 "Design of High Speed MOS Multiplier and Divider Using 
Redundant Binary Representation," IEEE Proceeding of 
Symposium on Computer Arithmetic, 1987, pp. 80-86). This 
algorithm considerably shortens the time for . remainder 
subtraction by using carry-propagation-free SD addition. 

15 However, it is much more complex because in each iteration 
the SD algorithm must check three most significant digit 
(MSD) bits of the remainder to decide the quotient digit in 
the set of {-1, 0, 1}, .and then perform the SD addition. 
Moreover, the final SD result must be converted to binary 

20 representation. Also note that the signed-digit addition is 
more complicated than the conventional carry-save adder 
(CSA). 

Another type of algorithms entirely avoids the slow 
25 subtract-detect-shif t type of operation previously 
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mentioned. They transform the division operation to a series 
of multiplication operations that converge to the original 
quotient. Among the examples are the constant convergence 
(S. Waser and M. J. Flynn, Introduction to Arithmetic for 
5 Digital Systems Designers, New York: CBS College Publishing, 
Chap. 5, 1982) and quadratic convergence (P. Markenstein, 
"Computation of Elementary Functions on the IBM RISC 
System/6000 Processor," IBM Journal of Research and 
Development, Vol. 34, 1990, pp. 111-119; D. A. Patterson and 

10 J. L. Hennessy, Computer: A Quantitative Approach, San 
Mateo, CA, Morgan Kaufman, 1990) division algorithms which 
are based on Newton-Raphson algorithm. They are often found 
in multiplier-based processors. They are still sequential 
type of operation to certain degree, and obviously require 

15 much more shif t-and-add operations. 

There is an on-line division algorithm that facilitates 
serial/serial division operation (K. S. Trivedi and M. D. 
Ercegovac, "On-Line Algorithms for Division and 

20 Multiplication," IEEE Trans, on Computers, Vol. C-26, No. 7, 
July 1977). This algorithm has advantages such as that: (a) 
it is pipelined at digit level; (b) all operands and results 
are communicated digit serial, and (c) result digits are on- 
line obtained after a few initial delay. On the other hand, 

2 5 among some of its disadvantages are: (a) it requires more 


3 


4 


4 

0438-0916-Q:SHC 


complex three-input signed-digit addition operation; (b) it 
needs more complicated quotient decision circuit for range 
detection of the remainder, and (c) output results have to 
be converted to binary representations. 

SUMMARY OF THE INVENTION 


10 In this work, a fast radix-2 division algorithm and its 

architecture is proposed* The algorithm adheres to the 
shif t/subtract-and-add type of division operation. The key 
idea behind this algorithm is to separate the, sign detection v 
operation of remainder from the remainder subtraction 

15 operation. By taking the absolute values of the remainders, 
we can successively subtract the remainders without the need 
of knowing the signs of remainders, while signs of the 
remainders can be decided in parallel and independently at 
the same time. To enhance the algorithm's performance, 

20 several design techniques were incorporated into its 
architecture realization. 


The new algorithm and its architecture - try - to - retain^, 

^ ^ the moot — advantages of the mentioned algorithms as possible, 

25 and simultaneously gets rid of their disadvantages. The 
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algorithm adopts non-restoring division operation and CSA 
type of operation for fast subtraction. Quotient digit set 
of {1, -1} is assumed for fast quotient conversion to binary 
representation. The algorithm is also an on-line algorithm 
that facilitates highly pipelined operation while it is much 
simpler than the existing on-line algorithms. 

This object of the present invention •,«HFej^ fulfilled by 
providing a method for finding a quotient Q = a^a^a^. . ... a^^ 
from a divisor Y = y,y, ....y and a dividend X = xx^....x. 
The method comprises the following steps of: (a) aligning 
the first non-zero bit of X with the first non-zero digit of 
Y; (b) defining a signed-digit partial remainder series R. 
where = Y, a first sign series of the partial remainder S. 
where S^^ = 0, a second sign series of the partial remainder 
S^., and a counter i beginning from zero; (c) subtracting X 
from R. which yields next signed-digit partial remainder 
R.^^; (d) setting the sign of R.^^ to S^.^^; (e) setting the 
result of exclusive-OR of S. and S^^, to the true sign of the 
next remainder S. , ; (f) setting a, to 1 if S.^, = 0 or 
R^^ = 0; (g) setting a. to 0 if S.^^ = 1; (h) inverting the 
signs of all digits of R^^ if S.^^ = 1; (i) shift R^^ left by 
one bit; (j) adding 1 to i ; and (k) repeating said steps (c) 
to (j) until i reaches a predetermined value or R.^^ = 0. 
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Further scope of applicability of the present invention 
will become apparent from the detailed description given 
hereinafter. However, it should be understood that the 
detailed description and specific examples, while indicating 
5 preferred embodiments of the invention, are given by way of 
illustration only, since various changes and modifications 
within the spirit and scope of the invention will become 
apparent to those skilled in the art from this detailed 
description . 

DETAILED DESCRIPTION OF THE EMB nnTMTCMTfi 

15 For achieving fast division, a new division algorithm 

is proposed and discussed hereinafter. 

NEW DIVISION ALGORITHM 

20 Given normalized n-bit sign magnitude operands 

1/2<|X|<1, 1/2<|y|<1 (this limitation is used to automatica 
lly align the first non-zero bit of X with the first non- 
zero digit of Y in the later discussion, which is avoided in 
circuit realization), quotient of Y/X can be solved using 

25 the following principles, where the quotient digits 
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Q2=a^aj^- a^a^. . .a^ is in sign-magnitude representation and a^ is 
the sign bit. 

Principle 1 : 

5 a^ equals to the result of exclusive-OR of the sign 

bits y and x of Y and X, respectively, i- e. , a = y. ^ x . 

Principle 2: 

Partial remainder R.^, can be solved by modifying the 
10 conventional, nonrestoring algorithm stated by K. Hwang as 
follows. The signed-binary quotient Q has its quotient digit 
g. € { 1 , -1 } , and 

R^, = 2|R - g^^/X| (1) 

15 

where R^= | Y | , ^0=1, g/ is the i-th pseudo quotient digit. 
Since R.^, is always positive, Eq. (1) can be rewritten as 


R^, = 2|R - X| 



f 1, if = 0 

1' if = 1 
-1, if S,^, = 1 


^i = 


L 0, if 


.f g. = -1 

^i + 1 


(2) 
- (3) 


(4) 


The sign of remainder (R. - X) , 
True sign of i-th remainder = S^^, ® S^^ , 
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Z. = Zero Flag, Z^^ = 0 , and 
So = ^ro = Sign{R^} = 0. 

The algorithm can be performed using the method 
5 described below • 

NEW DIVISION METHOD 

Define signed magnitude numbers = YiYj* 

10 = x^ XjX^ x^, and = q^q^, * q^q^ q^ in signed-binary 

representation, q^€{-l, 1}, and converted to sign magnitude 
representation = a^^o* ^1^2* • • - ^b' a. ^ {0, 1} . The quotient 
Q2 of xyy^ can be solved by the following steps: 

15 Step 1: 

a = y ® X • 

s -* s s 

step 2 : 

Define Y = y^y^ y^, X = x^x^ x^, Q = a^a^a^ a^, 

20 R^j = Y, i == 0, and S^^ = 0, 

Step 3: 

Subtract X from R. and yield next signed-digit partial 
remainder R. - Set the sign of R. , to S^. , (note that the 
25 sign of R.^^ equals the sign of first non-zero digit of R^j)- 
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Set the result of S. ® S^+i true sign of the next 

remainder S. Set a. to 1 if S. . = 0 (means the remainder is 
positive) or R.^^ = 0. Set a. to 0 if S.^^ = 1 (means the 
remainder is negative) • 

5 

Step 4: 

If S. =1, then take the absolute value of R. , (by 

i+i ' i+i ^ 

inverting the signs of all digits). Shift R.^^ left by one 
bit. -Add 1 to i. Repeat step 3 until i reaches a 
10 predetermined value or R^^ = 0. 

For better comprehension, two examples are used to 
demonstrate the division method hereinafter: 
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EXAMPLE 1: 


y = oioioooi^ = 81 


X = 00001001^ = 9 (X will be shifted left 3 places to 


align its highest non-zero digit with highest non-zero digit 


1010001 Rq = Y 

-1001000 X 


10 


Q> 0001001 > 0 then S^, = 0, 


S = S^, ® S = 0, a = 1 


SHIFT LEFT ONE BIT 

0010010 

15 - 1001000 


TAKE ABSOLUTE VALUE 
20 0110110 
SHIFT LEFT ONE BIT 

1101100 
- 1001000 X 


-0110110 < 0 then S^^ = 1 , 

S^ = S^2 ® Sj = 1, a^ = 0 


25 0100100 R3 > 0 then S^^ ^ 0 , 

S3 = S e S^ = 1 , a^ = 0 

SHIFT LEFT ONE BIT 

1001000 
- 1001000 X 


30 


0000000 R = 0 then a = 1 

4 3 

Result: 

The quotient = a^a^a^a^ = 1001^ = 9, 
and remainder = 0 


10 

// 
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EXAMPLE 2: 

In this example, note that a digit with a bar means a digit 
having negative value, e. g. 


11 = 2 + (-1) = 1. 


Y = lOlOlllOj = 174 
X = 11^ = 3 


10101110 = y 


-11 

10 


01101110 < 0 then S^^ = 1 , 

TAKE ABSOLUTE VALUE 

01101110 

15 SHIFT LEFT ONE BIT 

11011100 

- 11 X 


10011100 R^ < 0 then S^^ = 1, 
20 = S^^ e Sj = 0, a^ 

TAKE ABSOLUTE VALUE 

10011100 

SHIFT LEFT ONE BIT 

100111000 

25 . - 11 X 


01111000 R3 > 0 then S^3 = 0 , 

S, = S^ © S, = 0, a, = 1 

SHIFT LEFT ONE BIT 
30 11110000 

- 11 


00110000 R^ > 0 then S^^ = 0 , 

S^ =, © S3 = 0 , a3 - 1 

35 SHIFT LEFT ONE BIT 

01100000 
- 11 X 


10100000 R < 0 then S- = 1, 
40 S' = S^^ ® S^ = 1, a^ = 0 

TAKE ABSOLUTE VALUE _ 

10100000 


11 
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SHIFT LEFT ONE BIT _ 

101000000 

- 11 X 


00000000 R =0 then = 1 , = 0 

6 So 


The quotient = a^a^a^a^a^a^a^ = 0111010^ = 58, 
and remainder = 0 

10 

Since absolute values of the partial remainders are 
computed instead of their actual values, the algorithm 
facilitates parallel computations of partial remainder and 
quotient digit. To further speed up the operation of 
15 subtraction in the preferred embodiment, we use specified 
signed-digit operation. 

SPECIFIED SIGNED-DIGIT SUBTRACTION 

20 Since computations of Eq. (2) involves only the 

subtraction operation of two positive numbers, R. and X, we 
can speed up the computation by defining the CSA-like 
operation as follows. 


25 y . - X. = 2c.^^ + t. (5. a) 

t. - 

wherein y., r. ^ {-1, 0, 1} 


t. + c. = r. (5.b) 

J J J 
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X., t e {0, 1} 
c. e {-1, 0}. 

Here, a signed-digit y. (represents the j-th digit of 
5 R.) subtracts a binary digit x., then generates carry c.^^ and 
intermediate result t.. The finial result r. (represents the 
j-th digit of R.^^) is obtained by adding t. and the carry-in 
bit c.* Since r. e {-1, 0, 1}, there will be no carry 
generated from t. + c.. As a result, the specified 
10 signed-digit subtraction efficiently eliminates carry 
propagation* In addition, the complexity of this operation 
is similar to that of conventional CSA. Example 3 depicts 
the modified subtraction method where T. = tt, ....t and 

1 12 n 

C". "~ • ■ • • • 

I 12 n 


13 

14- 


0438-0916 


-Q:SHC 


EXAMPLE 3: 


10 


15 


20 


25 


30 


35 


40 


Y = OlOlOOOlj = 81 


X = 00001001, = 9 



1010001 
-1001000 

0011001 
-0001000 

0001001 


SHIFT LEFT ONE BIT 


0010010 
- 1001000 

1011010 
-1001000 

Tiooioio 


TAKE ABSOLUTE VALUE 


11001010 


SHIFT LEFT ONE BIT 


110010100 

- 1001000 

111011100 
-011011100 

001100100 

SHIFT LEFT ONE BIT _ 

11001000 

- 1001000 


10000000 
-01000000 


T. 


> 0 then S^^ = 0 , 


S = S 


n 


a. = 1 


< 0 then S^2 ^ ^ ' 
e Sj = 1, 


s = S^ 


= 0 


X 

T3 

R^ > 0 then S^3 = 0 , 

= © S, = 1, a, = 0 

X 


R, = 0 then a, = 1 


Result : 


The quotient = 1001 = 9, and remainder = 0 


14 

IS 

* 


ft 


0438-0916-Q:SHC 


As shown in the above example, T. and C. are calculated 
first, then R can be easily decided. The .truth table of t. 
and c^^^ value are listed in table 1, where the signed-digit 

represented by two bits , r. 


of 


R 


is 


and 


r.' = sign{r.). r.^ = 


TABLE 1: 


TRUTH TABLE OF t. AND c. 


X 

J 

r.^ 

r.^ 
J 


t. 

0 

0 

0 

0 

0 

0 

0 

1 

0 

1 

0 

1 

1 

1 

1 

1 

0 

0 

1 

1 

1 

0 

1 

0 

0 

1 

1 

1 

1 

0 


Note that r. represents j-th digit of R. here. 


The truth table of r. (represents the j-th digit of 


R ,) is listed in table 2 
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TABLE 2 


TRUTH TABLE OF r. (the j-th digit of R^^ ) 


c 
J 

t. 

J 

i 


0 

0 

0 

0 

0 

1 

0 

1 

1 

0 

1 

1 

1 

1 

0 

0 


15 


20 


From the above tables we can see that the signed-digit 
subtraction can be achieved by a simple digital circuit, and 
can be pipe-lined for better performance. 

CONCLUSION 


25 
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In summary, the division algorithm have the advantages 
as follows: 

a) It uses a smaller quotient digit set of {1, -1} than 
{-1, 0, 1}, that simplifies that quotient decision circuits 
like some known algorithms do, but achieves the exact 
division and trivial conversion of the results from signed- 
binary representation to binary representation. 
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b) It needs no quotient estimator. 

c) In each iteration, the algorithm computes partial 
remainders without knowing the signs of previous remainders 
and decides the signs of remainders independently and in 
parallel. In addition, these two operations are done in 
pipelined fashion and in digit level with maximum throughput 
rate . 

d) Its architecture *d-9- basically consists of the simple 
CSA type cells. 

e) It can handle either positive or negative operands. 

From above discussion, the proposed division algorithm 
and its architecture is very efficient. The new algorithm's 
realization is composed of highly regular cellular array, 
which is suitable for VLSI implementation and can be easily 
extended to bit-parallel implementation .i^C O 


The algorithm can be extended to higher radix divisions 
such as radix-4 division. Since the remainders are taken 
absolute values, the digit set contains only digits 1 and 2 
is sufficient for. the entire radix-4 operation. This greatly 
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reduces the number of search regions for the quotient 
digits , in contrast to the bigger set of { 0 , 1 , 2 , 3 } that 
existing algorithms allow, 

5 While the invention has been described by way of an 

example and in terms of several preferred embodiments, it is 
to be understood that the invention need not be limited to 
the disclosed embodiment. On the contrary, it is intended to 
cover various modifications and similar arrangements 
10 included within the spirit and scope of the appended claims, 
the scope of which should be accorded the broadest 
interpretation so as to encompass all such modifications 
and similar structures. 
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